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By critical analyses of the order parameter of symmetry breaking, we have researched the phase 
transitions at high density in D = 2 and D — 3 Gross-Neveu (GN) model and shown that the 
gap equation obeyed by the dynamical fermion mass has the same effectivenesss as the effective 
potentials for such analyses of all the second order and some special first order phase transitions. 
In the meantime we also further ironed out a theoretical divergence and proven that in D = 3 GN 
model a first order phase transition does occur in the case of zero temperature and finite chemical 
potential. 
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The symmetry restoring phase transition at high temperature and high density has been an interesting research 
topic which is closely relative to evolution of early universe and phase transition of the nuclear matter [§ §• 

A good laboratory to research the phase transition is the models of four-fermion interactions including the Gross- 
Neveu(GN) models in D = 2 and fl = 31 time -space dimensions. On the phase transitions in these models, extensive 
investigations have been made || [?], ||, ^ [H^, 11, |l2|, [l3|, some common conclusions have been reached but there is 



X 



still particular question, e.g. whether first order phase transition exists or not in D — 3 GN model, which needs to 
be further clarified. 

For determining the order of a phase transition, one can analyze critical behaviors of dynamical fermion mass as the 
order parameter of symmetry breaking at a critical temperature and/or at a critical chemical potential. The phase 
transition will be second or first order if the dynamical fermion mass goes to zero at the critical point continuously 
or discontinuously. 

The fundamental approach of the critical analyses should be using an effective potential. However, in many cases, 
the gap equation is a more simple and direct means. So called the gap equation we mean the essential part of the 
Schwinger-Dyson equation obeyed by the dynamical fermion mass as the order parameter of symmetry breaking. It 
can be independently derived from the fermion condensates induced by the four-fermion interactions without using an 
effective potential, though it has been proven that the gap equation can also come from the extreme value condition 
of an effective potential [ jO] . If a phase transition is second order and this normally happens when the gap equation 
has the only non-zero solution, then the critical analysis can be made simply by means of gap equation. This has been 
carried out in the demonstrations on high temperature second order phase transition for fixed chemical potential in 
D = 4 Nambu-Jona-Lasinio model JT^] ,Z) = 3 JL6) and D — 2 JTt) GN model. However, if a phase transition is first 
order and this is usually the case when the gap equation has two or more non-zero solutions, then how should one 
make the critical analysis ? Whether or in what conditions can one do that simply based on the gap equation ? 

To answer the above questions, in this paper we will generalize critical analyses to the case including first order 
phase transitions. As far as the four-fermion models are concerned, it is most possible that the first order phase 
transition occurs in the case of high density. Hence we will investigate symmetry restoring phase transitions at high 
density for a given temperature, especially for zero temperature case, in D = 2 and D = 3 GN models. We will use 
both gap equations and effective potentials, and specially examine effectiveness and limitations of the gap equation 
for critical analyses of first order phase transition. 
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The Lagrangian of the models will be expressed by 



N N 
k=l k=l 



where ipk{x) are the fermion fields with N "color" components and g is the four-fermion coupling constant. The 
discussions will be made in the fermion bubble diagram approximation which is equivalent to the leading order of 
the 1/N expansion. It is assumed that higher-order corrections in 1/N do not change the qualitative feature of the 
models, since our models involve only some discrete symmetries which will be spontaneously broken or restored [ fl8[ . 
Our basic starting point is the gap equations coming from the fermionic bilinear condensates induced by the four- 
fermion interactions in Eq.(l.l), besides the corresponding effective potentials. The gap equations have been given 
in Ref. ]l6|, |l7| . The corresponding effective potentials may be obtained from the relations between the gap equations 
and the effective potentials proven in Ref. pA L Here we will directly use the known explicit results of them. 

The paper is arranged as follows. In Sect. □] we discuss D = 2 GN model, including analyzing critically the second 
order high density phase transition at high temperature T by the gap equation and in particular, the first order high 
density phase transition at T = by the gap equation combined with the corresponding effective potential. In Sect. 
Ill, similar discussions to ones in Sect, [n] will be made in D = 3 GN model. The first order high density phase 



transition at T = in this model will be proven respectively by critical analyses based on the gap equation and the 
effective potential. A difference between it and general first order phase transition and this difference's significance 



for usefulness of the gap equation analysis will be indicated. Finally in Sect. IV we come our conclusions 



II. D = 2 GROSS-NEVEU MODEL 



In this model [p"7j], there exists a matrix 75 = — a 3 anticommunicating with 7 = a 1 and 7 1 = ia 2 , where 
— 1,2,3) are the Pauli matrices. By the Mermin- Wagner- Coleman theorem, a continuous symmetry cannot 
be spontaneously broken in 2 dimensions |19| . Dynamical generation of the fermion mass will spontaneously break 



only the discrete chiral symmetry \d ■ ip(t,x) 75^(4, a;) and the special parity V\ : tp(t,x) — l -> ^^(t 
Hence possible phase transition at high temperature and high density will relate to restoration of only the \d and 
V\. However, in 2 dimensions, symmetry restoration could occur only in the mean-field approximation |3| ^| which 
is equivalent to fermion bubble diagram approximation [T^]. Because the former can be used only if N — ► 00 or N 
is large and only a finite segment of one-space dimension system is considered [^0|, here it is assumed that the use 
of the fermion bubble diagram approximation will be in the above conditions. In 2 dimensions, Lagrangian (1.1) 
is perturbatively renormalizable. Assume the four-fermion interactions in Eq.(l.l) could lead to the condensates 

Sfc=i(V' fc ' ! / ; fc) 7^ when T = fj, = and the thermal condensates J2k=i (r ^&)t 7^ when T ^ 0, where /i is the 
chemical potential of fermions, then we will obtain the renormalized gap equations in respective cases and finally 
reduce the gap equation at T 7^ to the following form 

raff)) 

In = F 1 (T, fi, m) = h(y, -r) + h (y, r) (2.1) 

777 

with the denotations 

I 1 (y, T r)= — , (2.2) 

Jo v 2 - + V cxp(^x z + y z =F r) + 1 

where y — m/T, r — n/T, m(0) and m = m(T, /_i) are respectively the dynamical fermion masses at T = fj, = and 
at T ^ and/or // 7^ 0. It is easy to verify that Fi(T, fi,m) is an increasing function of T and /j, but a decreasing 
function of 777, and from Eq.(2.1), we can assume that at some critical temperature T c and critical chemical potential 
/i c , the dynamical fermion mass to becomes zero. This will lead to a critical equation. At a high temperature T, by 
using the high temperature expansion of F\(T, /i, to) Jl5[ we may reduce the gap equation (2.1) to 

ln ^ + 7 = 9(T, + h(T, + O [(£) 4 ] , 

g(T, fj,) = |C(3) - 93C(5)(£) 2 + ifi C (7)(^) 4 , (2 . 3) 
h(T, fj,) = 7C(3)(£) 2 - 31C(5)(£) 4 + 127C(7)(£) 6 , 



-x 
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where C( s )( s = 3,5,7) are the Riemann zeta functions and 7 is the Euler constant. Noting that the lnm term in 
Eq.(2.1) has been cancelled, thus we can set y = in Eq.(2.3) and obtain the critical equation 

m ?Z^) +7=/l( T c , Mc ). (2.4) 

In the derivation of Eq. (2.4) an implicit fact is that Eq. (2.3) has the only solution m and m can go to zero 
continuously as T and [i vary. This is just the characteristic of a second order phase transition. In view of this, Eq. 
(2.4) can be called the critical equation of a second order phase transition. The critical points determined by Eq. 
(2.4) are all second order ones. For verifying this conclusion, we can discuss from Eqs.(2.3) and (2.4) critical behaviors 
of the squared order parameter m 2 at \x c when T is fixed. Replacing T c by T in Eq.(2.4) and substituting it into Eq. 
(2.3), we get 

™ 2 - (Mc ™ fi?)f(T,Vc)/g(T,iJ,), when^~/x c , (2.5) 

where 

/(T> c ) = 7C(3) - 62C(5)(£^) 2 + 381C(7)(£|0 4 . (2.6) 

Eq.(2.5) shows that the squared order parameter m 2 goes to zero continuously when /i — > /i c in the form of [i c — /x, 
hence for a given high temperature T, the symmetry restoration phase transition at [i c is second order indeed. 

Eq.(2.4) comes from high temperature expansion of F±(T, fJ.,m) and is inapplicable to discussion of high density 
phase transition at low temperature and T = 0. For discussion of the latter, we have to use the effective potential. 
In the following we will only focus on the high density phase transition at T = 0. Let us first go back to Eq.(2.1). 
Through changing the integral variable by z = (x 2 /y 2 + l) 1 / 2 we can express 



Fi(T,n,m) = J 



dz 



where a = — . In the limit T — ► or y — ► 00, we have 



1 

+ (—a — > a) 



e y(z-a) + I 



(2.7) 



Fi(T = 0,n,m) = e(n-m)ln ^ + ^ — . (2.8) 

m 



Substituting Eq.(2.8) into Eq.(2.1) we obtain 



which has the solutions 



and 



, m(0) . , u + \l u 2 — m 2 

In — ^ = - m) In ^- (2.9) 



m = m(0), when /j, < m (2-10) 



{to(0)[2^-to(0)]} 1/2 = mi, 
m(0) 

when /i > m and — - — < fi < m(0). (2-11) 

,(2) 



On the other hand, the effective potential V e * i (T, fx, m) has been proven to satisfy the equation |14| 



dm 7r 



In— 77^- + Fi(T,^,m) 
m[0) 



(2.12) 



(2) 

noting that the gap equation (2.1) multiplied by m is just the extreme conditions given by dV£> f {T,n,m)/dm = 0. 

(2) 

Hence, when T = 0, V^J^iT = 0, /i, m) for m(0)/2 < fx < m(0) will have the three possible extreme points m = 0, mi 

and m(0). To determine which one of them corresponds to a minimum or a maximum of V^y^ (T = 0, /i, m), we may 
calculate and obtain directly from Eq. (2.12) 

d 2 V^ f (T = 0, fi, m) _ 1 j ln 4w + 1 when M ^ m 

0^2 -~ \ ln t±Vj£_^ aL- when M >m' (2 ' 13) 

^ m (0) ( A1 +^ /i 2_ m 2)^ M 2_ m 2 ^ 
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where Eq. (2.8) has been used. Eq. (2.13), combined with Eqs. (2.10) and (2.11), give that 

d 2 vg}(T = 0,^m) 



m—O 



7r ml 



< when fj, < m(0)/2 
(0) \ > when fx > m(0)/2 ' 



(2.14) 



d 2 V$(T = 0,Li,m) 



dm 2 



1 2 M -m(0) m(0) 

< 0, when — — < fj, < m(0 

7r m(0) — /x 2 



(2.15) 



and 



c> 2 y e ( g(r = o,M,m) 



9? 



> 0, when fi < m(0). 



,(0) 



(2.16) 



Eqs. (2.14)-(2.16) indicate that 

1) When [i < , V^yl(T = 0,fi,m) will have the maximum point m = and the minimum point m = m(0), thus 
the model has the same spontaneous symmetry breaking as one at /i = 0. 

2) When < M < m (0), K//( T = °> 

/x,m) will have three extreme points: m = (minimum),™ = mi (maximum) 
and m = m(0) (minimum), and this is just the characteristic of possible emergence of a first order phase transition. 

It is emphasized that in the above analyses we only use the known relation (2.12) between the gap equation and 
the effective potential and the explicit expression of the gap equation (2.1), and did not use an explicit expression of 

(2) (2) 

\Qff(T, /i, m). However, for determining the critical point, we have to find out the explicit expression of \Q^i(T, /i, m) 
from Eq.(2.12). The result is Q 



7Tb 

0,M,m) = — 



In 



m(0) 



2, M+ Vm 5 
m in 



EL 

2tt 



mf 



±9(» -m) 



(2.17) 



Assuming that m ^°- ) < ii < m(0), then for m = m(0) among the three extreme points, we will get 



y e ( / 2 }[T = 0, Ai ,m = m(0)] = -^[2 M 



47T 



2 m 2 ( 



<0, 



Hence, if /i < m(0)/y/2, m 

( 2) 

minimum V^y^(T = 0, fi,rr, 
transition point should be determined by 



when /i < — m(0). 
V 2 

m(0) will be the global minimum point of V^yi(T 



(2.18) 



and V^T = 



0, /x, m 



0,/i, m), because in this case, the 
mi) is a maximum. Obviously, the first order phase 



K// [T = 0, M , m = m(0)] = K// (T = 0, /z, m = 0) = 



— t/(2) rT — 



(2.19) 



and this will lead to 



11 = V2 m ^ = 

r(2). 



(2.20) 



i.e. when /i crosses over /x c , the global minimum point of \Qf {(T = 0, fx, m) will jump from m(0) down to through a 
barrier tunneling between the minimum points m = m(0) and m — 0. The critical chemical potential fi c of the first 
order phase transition at T = in D = 2 GN model in the leading order of l/N expansion given by Eq. (2.20) is 

consistent with the one obtained by other authors @, ||, [| [n| [ll} [l^| . Since the global minimum of V^iT = 0, /i, m) 
is also at m = m(0) when /i < m(0)/2, we can write the critical behavior of the order parameter m by 




when 
when 



^ < ^m(0) 
H > ^m(0) 



(2.21) 
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The above discussions show that for critical analyses of the order parameter in a general first order phase transition, 
it is not sufficient only using the gap equation, we must use the gap equation combined with corresponding effective 
potential. This is also applicable to critical analyses of the phase transitions at low temperature T 7^ in D = 2 GN 
model which are known to be first order. The effective potential at T 7^ may be derived from Eq.(2.12). By means 
of it and the gap equation (2.1), we will be able to obtain other results of the phase transition at low temperature in 
this model including location of the tricritical point etc., following similar methods to ones used in Ref |(| |7[ [L2| . 



III. D = 3 GROSS-NEVEU MODEL 



In this case Jl6| , the dynamical fermion mass will spontaneously break the special parities V\ and V 2 . ip{t t x 1 , x 2 ) 

T 



— l -y y 1 ^, —x l , x 2 ) and ip(t, x 1 , x 2 ) j 2 ip{t, x 1 , —x 2 ) and the time reversal T: ip(t, x 1 , x 2 ) 
After the gap equation at T = fi = is substituted, the gap equation at T 7^ becomes 

m(0) = m + F 2 (T, /1, to) 

with 

F 2 (T,/i,m) = T {in [l + e -( m ~' 1 )/ T ] + (-/i -> M )} , 

and the critical equation of second order phase transition will be 

to(0) = T c [ln(l + e^/ T -) + m (l + e -^ /T ") 

which determines critical points of second order phase transitions. For example, for a given finite T, when \x 
to w 0, so we will have m/T < 1. Then Eq.(3.1) with (3.2) will lead to 



7 2 V(-i,x 1 , a; 2 )|l). 

(3.1) 

(3.2) 
(3.3) 



m(0) = rin(l + e^ /T ) +Tln(l 
m 

o 



1 



TO 

T3 



2T[l + cosh(^/T)] 
Substituting Eq.(3.3) with T c replaced by T into Eq.(3.4) we get 

m 2 = 2Tsinh(/Lt/T)(/x c - n), for T ^ 



(3.4) 



(3.5) 



Eq.(3.5) shows that for a given finite temperature, the symmetry restoration phase transition at [i c is second order. 
The same has been shown for the phase transition at T c for a given /1 | ft6| . These conclusions are valid for T 7^ 0. For 
T = 0, we must go back the gap equation (3.1). We note that when fj, < m, lim F 2 (T, /i, to) = and by Eq.(3.1) this 

implies that m 



lim F2CT, a, to 

T->0 



to(0), if fi < m, or furthermore, to = to(0), if /J, < to(0). On the other hand, if > m, then 
/i — to. Substituting this result into Eq.(3.1) we will obtain that to(0) = n, but the presupposition 

leading to this result in fact should be fj, > m = m(0), the both apparently contradict each other. This indicates that 
when fi > m = to(0), the gap equation (3.1) can not be satisfied, thus there is no longer symmetry breaking and the 
order parameter to must be zero. The above results can be summarized as 



TO 



to(0) when fi < to(0) 
when fi > to(0) 



(3.6) 



Furthermore, when T 



0, the derivative dm/dfj, will be 
dm 



lim 



lim 



m— m(0) 

dF 2 (T,n,m) 
9/i 



dF 2 (T,n,m) 



■ lim e -M0)-rt/T 

T-»0 



^ m J m=m(0) 

0, when < to(0) 

— 1, when = to(0) 

-co, when > to(0) 



(3.7) 
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(noting that when fi > m(0), the gap equation (3.1) is no longer valid). Eqs.(3.6) and (3.7) show that when fx < m(0), 
m does not change as fx varies and keeps its value m(0) at fx = 0; however as soon as fx crosses over m(0), dm/dfi 
will change from into — oo and m will jump discontinuously from m(0) to zero. This implies that when T = the 
phase transition at the critical chemical potential fx c = m(0) is first order, i.e. (T, /i) = (0,to(0)) will be a first order 
phase transition point in T — fx phase diagram. On the other hand, we note that (T, fx) — (0, m(0)) is a solution of 
the critical equation (3.3) corresponding to second order phase transition when we take the limit T — > in Eq.(3.3). 
Hence we may reasonably conclude that (T,fx) — (0,m(0)) will be a tricritical point in the T — fx phase giagram of 
D = 3 GN model. 

This conclusion coincides with the one in Ref. [jl3| but different from the one in Ref. |Q where the authors claimed 
that there is no first order phase transition in D = 3 GN model. To clarify this problem, we can make further analysis 
of the phase transition at T = in the model by means of the effective potential which may be derived from the gap 
equation, as was made in D = 2 GN model, ft has been proven that the extreme value condition of the effective 
potential V^](T,fx,m) is @ 

eff \ =0 (3.8) 

dm 

with 

dv£){T,fx,m) m 

dm =^[m-m(0)+F 2 (T, fx, m)} , (3.9) 

noting that Eq.(3.8) is just the gap equation (3.1) multiplied by m. From Eq. (3.9) we may obtain V^)(T,fx,m) 
satisfying the condition V^^(T, fi, m = 0) = expressed by 

V e { V(T = 0, M , m) = i- j^Vox - m) - m(0)] + 0(m - M ) + . (3.10) 

For fixed fx, from Eq.(3.10) or the T -> limit of Eq.(3.9),we find out 



dV e f f (T = 0^,m,) 



in 



and 



dm 2tt 

d 2 V e f f (T = 0^,m) _ i 
9m 2 2ir 



= —{[/*- m(0)]6>(^ - m) + [m- m(O)]0(m - /x)} (3.11) 



{[/x - m(O)]0(/z - m) + [2m - m(O)]0(m - fx)}. (3.12) 



Variation of V^f|(T = 0, fx, m) as fi increases may be discussed by Eqs. (3.10)-(3.12). 

(3) 

1) fx < m(0). In this case, V^^[T — 0, fx < m(0), m] will have the maximum point m = 0(m < fx) and the minimum 
point m = m(0)(m > fx). It is indicated that location of the minimum point keeps at m — m(0) which does not 
change as fx increases. However, the minimum of the effective potential at m — m(0) 

V e f f [T = 0,fx< m(0),m]\ m=m(0) = ^-[fx 3 ~ m 3 (0)] (3.13) 

will rise as fx increases. Nevertheless, these results only show the same spontaneous symmetry breaking as one at 
T = fx = 0. 

2) fx = m(0). In this case, the effective potential will take the following form 

V} 3 1 [T = Q,fx= m(0), m] = ^-9[m ~ m(0)] 
1 3 Ztt 



to 3 (0) m 



2 2 m(0) 



(3.14) 



which indicates that in the total real axis segment < m < m(0), V^yifT = 0,/^ = m(0),m] = 0, thus the total 

(3) 

segment < m < m(0) including the original m = m(0) when ^ < m(0), will be the minimum points of y e ^[T = 
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0, /i = m(0),m]. This situation corresponds to that when fj, = m(0), the T — > limit of the gap equation (3.1) has 
infinite non-zero solutions. 

("3") 

3) (U > ra(0). In this case, dV^^[T = 0,/x > m(0), m}/dm — will have the only solution m = and by Eq. (3.12), 



d 2 V^}[T = 0,//> m(0),m] 



dm 2 



m=0 



= — [/x-m(0)] >0. 



(3.15) 



Thus m = will be the only minimum point of V^y^ [T = 0, \i > m(0), m] and this implies that the broken symmetries 
will be restored. The critical chemical potential /i c should be determined by the condition 



d 2 V e f f (T. 



0,fi,m) 



dm 2 



= 0. 



(3.16) 



m=0 



i.e. at /i c , m = must change from being a maximum point into being a minimum point. The result is that fi c = m(0). 
The above discussions show that as /i increases and finally crosses over \i c — m(0), the global minimum point of 

the effective potential V^j(T = 0,/i,m) will jump from m(0) down to thus Eq. (3.6) is reproduced. The phase 
transition is first order indeed. The analysis based on effective potential leads to the same conclusion as the one only 
based on the gap equation, however the latter is obviously more simple and direct in this model. We also indicate 
that the first order phase transition at T = in D = 3 GN model has a distinct feature, i.e. the transit from broken 
phase to symmetry phase does not undergo a barrier tunneling, since when [i is increasing, as we have seen above, 
no barrier emerges from between the minimum points m — and m = m (0) of V£' f (T = 0, /z, m). This is apparently 
different from general first order phase transition, e.g. the one in D = 2 GN model discussed in Sect. [n| It is assumed 
that it is just this feature that makes one be able to determine the critical behavior of the order parameter in the first 
order phase transition in D = 3 GN model merely by the gap equation approach. 



IV. CONCLUSIONS 



In this paper, we have discussed symmetry restoring phase transitions at high density in D — 2 and D = 3 GN model 
by means of analyses of critical behaviors of the dynamical fermion mass as order parameter of symmetry breaking. 
We have used both gap equation and effective potential approach. In these approaches, the gap equations multiplied 
by the dynamical fermion mass are simply the Schwinger-Dyson equations obeyed by the dynamical fermion mass 
and the effective potentials will be derived from the gap equations based on the fact that a gap equation multiplied 
by the dynamical fermion mass can also come from the extreme value condition of corresponding effective potential. 
This shows that the Schwinger-Dyson equations obeyed by the dynamical fermion mass in fact contain the essential 
elements of the phase structure of the discussed models. We have found that for second order phase transitions and 
some specific first order phase transitions e.g. the one at high density and T = in D = 3 GN model in which the 
jumping of the order parameter does not correspond to a barrier tunneling in an effective potential, the gap equation 
analysis alone has the same effect as the effective potential one. However, for general first order phase transitions, 
e.g. the one at high density and T = and low T ^ in D = 2 GN model, we must combine the gap equation with 
the corresponding effective potential so as to obtain correct conclusions. 

Our analyses also reproduce some of known essential results of the phase transitions in D = 2 and D = 3 GN 
model. They include that in D = 2 GN model, the phase transitions at high density will be second order when T is 
high and first order when T —> 0; in D = 3 GN model, the phase transition at high density will be second order when 
T is finite and first order when T = thus one can conclude that (T, /i) = (0, m(0)) is a tricritical point. The latter 
further clarifies the theoretical divergence of that if first order phase transition exists in D = 3 GN model. 
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